Non-trivial spin-texture of the coaxial Dirac cones on the surface of topological 

crystalline insulator SnTe 
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We present first principles calculations of the nontrivial surface states and their spin-textures in 
the topological crystalline insulator SnTe. The surface state dispersion on the [001] surface exhibits 
four Dirac-cones centered along the intersection of the mirror plane and the surface plane. We 
propose a simple model of two interacting coaxial Dirac cones to describe both the surface state 
dispersion and the associated spin-texture. While the out-of-the-plane spin polarization is zero due 
to the crystalline and time-reversal symmetries, the in-plane spin texture shows helicity with some 
distortion due to the interaction of the two coaxial Dirac cones, indicating a nontrivial mirror Chern 
number of -2, distinct from the value of -1 in Z2 topological insulator such as Bi/Sb alloys or Bi2Se3. 
The surface state dispersion and its spin-texture would provide an experimentally accessible way to 
determine the nontrivial mirror Chern number. 

PACS numbers: 73.20.r, 73.43. Cd 



I. INTRODUCTION 

Since the discovery of time-reversal symmetry pro- 
tected topological quantum states in two-dimensional 
(2D) Hg(Cd)Te-based quantum well structures, and sub- 
sequently that of the Z2 three-dimensional (3D) topolog- 
ical insulators^ ^, an enormous effort has been dedicated 
to finding other novel materials, which could support 
non-trivial topological states. One particularly fruitful 
recent direction has been to explore quantum states in 
condensed matter systems, which are protected by the 
symmetries of the lattice, leading to the so-called topo- 
logical crystalline insulators (TCIs).^ A practical real- 
ization of a TCI has been the prediction of SnTe with 
an ideal rocksalt structure in which the mirror symme- 
try of the lattice ensures the existence of robust metallic 
edge states.^. This theoretical prediction was verified 
quickly via angle resolve photoemission experiments on 
Pbi_^Sn^Tei^^i and Pbi_^Sn^Sei^. Recall that a Z2 
topological insulator such as Bi2Se3/Te3, which is pro- 
tected by time-reversal symmetry, contains a single Dirac 
cone at the center of the [111] surface plane, exhibiting 
a linear dispersion and chiral spin texture predicted the- 
oretically as well as observed experimentall}'^. In con- 
trast to the 3D Z2 topological insulators, however, a TCI 
protected by mirror symmetry contains an even number 
of Dirac cones on crystal surfaces symmetric about the 
[110] mirror planes, and its topologically non-trivial state 
is characterized by a nonzero mirror Chern numbei!^^. 
The interest in understanding the bizarre surface states 



and their spin-textures in SnTe, the pristine phase of an 
archetype TCI, is thus obvious, both in its own right and 
as a way of gaining insight into the properties of related 
substitutional compounds. 

In this work, we report first-principles calculations to 
investigate the [001] surface states of SnTe in rocksalt 
structure. By including the spin-orbit coupling, two gap- 
less Dirac cones centered along f X (the other two lie 
along TY), a mirror line of the crystal symmetry, are 
found in the surface spectrum, a result of the nontrivial 
band topology due to crystal symmetry. We examine the 
complicated surface band characters as well as the non- 
trivial spin textures, and propose a simple model which 
consists of two interacting coaxial Dirac cones centered at 
X. When the interaction between the two coaxial Dirac 
cones is turned on, the two surface bands avoid each other 
and a gap opens up except on the mirror line, forming the 
gapless Dirac cones centered along TX. While the out-of- 
the-plane spin polarization vanishs due to C/^y symmetry 
as well as the superposition of the mirror and time rever- 
sal symmetries, the in-plane spin textures show helicity 
with some distortion due the interaction of the two Dirac 
cones. The overall spin texture reveals the nontrivial mir- 
ror Chern number with the value of -2, which is distinct 
from that of -1 i n Zg to pological insulators such as Bi/Sb 
alloy or Bi2Se3.'^'^3li4l ^]^g surface state dispersion and 
the associated spin-texture would provide an experimen- 
tally accessible way to determine the nontrivial mirror 
Chern number. 
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FIG. 1: (Color online) Crystal structure and surface 

bands of SnTe. a. Rocksalt crystal structure of SnTe. A 
[001] surface plane is shown, b. FCC Brillouin zone (BZ) 
of SnTe and the [001] surface BZ. The [1 -1 0] mirror plane 
crosses the 2D surface BZ along FX. c. Electronic structure 
of SnTe [001] surface around X. The surface bands are in- 
dicated by thick green lines and bulk bands by the shaded 
purple area. A Dirac point of surface state appears along the 
f X direction, while the surface band is gapped along XM. 



II. BAND STRUCTURE CALCULATIONS 

We first extract both the electronic band structure and 
the spin texture of the SnTe surface states from first- 
principles calculations, which was carried out within the 
framework of the density functional theory (DFT) using 
pseudo-potential projected augmented wave methocP^ 
as implemented in the VASP package^^. The general- 
ized gradient approximation (GGA)^^ was used to model 
exchange-correlation effects. The spin orbital coupling 
(SOC) is included in the self-consistent cycles. The sur- 
faces are modeled by periodically repeated slabs of 33- 
atomic-layer thickness with 13 angstrom wide vacuum 
regions, using a 12x12x1 Monkhorst-Pack k-point mesh 
over the Brillouin zone (BZ) with 208 eV cutoff energy 
The room temperature crystal structure of SnTe in an 
ideal sodium chloride structure was used to construct the 
slab without any rhombohedral distortion. The experi- 
mental lattice constant of SnTe with the value of 6.327 
A was used^^. The self-consistent Bloch wavefunctions 
associated with the surface states were decomposed into 
cubic spherical harmonic orbitals and projected on to 
various atomic sites. For each atomic site in SnTe, we 
obtained the charge density and the three components of 
the spin direction from the expectation values of the x, 
y and z Pauli matrices. 

In order to clarify the relevant structural aspects, SnTe 
in an ideal sodium chloride FCC crystal structure is 
shown in Fig. la. Fig. lb shows the FCC bulk Brillouin 
zone (BZ) as well as the [001] surface BZ. The high sym- 
metry points r, L and X in the bulk BZ are projected 
on the f , X and M in the surface BZ, respectively. The 



[1-10] mirror plane is perpendicular to the [001] surface 
plane and projected onto the surface BZ along the f X 
direction. The projected bulk band structure along the 
high symmetry lines f -X-M is shown in Fig. Ic by purple 
area. The surface bands are highlighted by thick green 
lines, isolated from other bulk bands. A gapless Dirac 
cone with the Dirac point sits along f -X on the two sides 
of the X point at the Fermi level, while the surface states 
along X-M are gapped. 

The complicated surface states of SnTe can be under- 
stood by a simpler picture, which consists of two Dirac 
cones. Fig. 2a shows two coaxial Dirac cones centered 
at the X point. The X and Y points, i.e. (0.5, 0) and 
(0, 0.5), are equivalent and we concentrate on the one 
at (0.5,0) hereafter. These two Dirac cones carry differ- 
ent band characters and can be distinguished through an 
analysis of the wavefunctions of the Sn and Te atoms on 
the surface layers of the [001] slab. The nature of the 
differences in potentials of the Sn and Te atoms sepa- 
rates the two Dirac cones in energy vertically as shown 
in Fig. 2a. If there were no interactions between the two 
Dirac cone states, the two Dirac cones will intersect each 
other in the 2D surface BZ in an ellipse, which only in- 
tersects the mirror line along f-X (blue line) with two 
points on the two sides of the X point. In Fig. 2b, when 
the interaction between the two Dirac cones is turned on, 
a gap opens up along the elliptical overlap region of the 
two Dirac cones, except at the two points protected by 
the mirror symmetry. 

In order to trace the band character of the SnTe sur- 
face states, we decompose their charge distribution into 
Sn and Te partial contributions presented in Figs. 2c-f. 
For the surface valence states (Fig. 2c), we find that the 
Sn fraction is dominant within an elliptical region be- 
tween the two Dirac points and centered at the X point, 
while in Fig. 2d the Te fraction is larger outside this re- 
gion. But, for the conduction surface states (Figs. 2e 
and f), this trend between the partial contributions from 
Se and Te atoms is reversed. To enhance visual clar- 
ity, we patch the charge fraction map onto the energy 
dispersion in Fig. 2g. These results of first-principles 
calculations qualitatively agree with the picture of two 
interacting Dirac cones shown in Fig. 2b and allow us 
to attribute the origin of the two coaxial Dirac cones to 
distinct atomic orbitals. One of the cones with opening 
to the high energy side is more Sn-like with p^-orbital, 
while the other cone with opening to the low energy side 
is more Te-like with p^^-orbital. The Sn-like cone has the 
X-point Dirac point lower in energy than that of the Te- 
like cone. 

In Fig. 3, we plot surface band energies using a color 
scale and the associated constant energy contours, to- 
gether with the in-plane spin-texture. The spin direction 
of each state was obtained here by calculating the ex- 
pectation of the Pauli matrices from the first 6 atomic 
layers from the top surface of the SnTe slab. For the 
valence surface states shown in Fig. 3b, the highest en- 
ergy state is located at the Dirac points lying along the 
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FIG. 2: (Color online) Two interacting coaxial Dirac cones. Schematic diagram for a non-interacting (interacting) coaxial 
Dirac cone centered at the X point is shown in a (b). As the interaction between the two Dirac cones is turned on, gaps open 
up except at the two new Dirac points on the two sides of the X point. The new gapless Dirac cones are protected by the 
presence of the mirror plane which intersects the surface BZ along FX (bold blue line). Maps of the fraction of the partial 
charges on Sn and Te atoms for the valence surface bands are shown in c and d, respectively, while those for the conduction 
surface bands are shown in e and f. Dashed lines are the constant energy contours, g. Surface band dispersion with colors 
indicating the fraction of partial charges on the Sn and Te atoms. 
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FIG. 3: (Color online) Non-trivial spin-texture of SnTe 

surface states. Spin textures of the conduction and valence 
surface bands are shown in a and b, respectively. Energy is 
represented by the color scale, and constant energy contours 
are shown. Green line indicates the FX direction and green 
dots denote the Dirac points. Schematic band structures for 
the conduction and valence surface bands along FX direction 
are shown in c and d, respectively. Solid lines are Te-like 
Dirac cones and the dashed lines are Sn-like. Red (blue) ar- 
rows indicate the spin along positive (negative) y-direction. 



FX mirror line. The constant energy contours close to the 
Dirac points are seen not to be of a perfect circular shape, 
which implies an anisotropic Fermi velocity around the 
Dirac points. Between two Dirac points, an energy valley 
is centered at the X point. The energy contours exhibit a 



Lifshitz transition^. As we go to energies below the Fermi 
energy, the constant energy surface changes its topology. 
The two disconnected hole pockets next to the X point 
at high energy become a large hole and a small electron 
pocket both centered at the X point at low energy. In 
the two coaxial Dirac cones, the large hole pocket is as- 
sociated with the Te-like lower Dirac cone and the small 
electron pocket is associated with the Sn-like upper Dirac 
cone. A similar change in the Fermi surface topology oc- 
curs in the conduction surface bands (Fig. 3a). At high 
energies above the interaction region, the large electron 
and small hole pocket centered at the X point are asso- 
ciated with the Sn-like upper Dirac cone and the Te-like 
lower Dirac cone, respectively. 

It is interesting to see how the spin-textures within 
our picture of two co-axial Dirac cones play out based on 
the spin-textures derived from our first-principles com- 
putations. Note that for the conduction bands shown in 
Fig. 3a, the spin texture shows counter-clockwise rota- 
tion around the X point and clockwise rotation far away 
from the X point. As shown in the schematic diagram 
of Fig. 3c, states near the X point are associated with 
the Te-like lower Dirac cone while those far away from 
the X point belong to the Sn-like upper Dirac cone. The 
two coaxial cones therefore should mimic similar chirality. 
They should both have clockwise spin rotation in the up- 
per cone and counter-clockwise in the lower cone. Along 
the f X direction, these two opposite spin states meet at 
the Dirac points. Along the XM direction, the spin polar- 
ization diminishes and switches direction around the Lif- 
shitz transition. The spin texture of the valence surface 
states can be understood along similar lines. A clockwise 
rotation in the inner region around X is associated with 
the Sn-like upper Dirac cone, while a counter-clockwise 
rotation in the outer region far away from the X point 



is associated with the Te-hke lower Dirac cone. We note 
that the spin chirahty of the surface Dirac cone in a typ- 
ical strong topological insulator like Bi2Se^is the same 
as that of the Sn-like and Te-like Dirac cones we have 
found here on the surface of SnTe, and bears a nontrivial 
mirror Chern number further discussed below. 



III. SIMPLIFIED k p MODEL 

In order to make our proposed two-coaxial-cone pic- 
ture more concrete, we now discuss a 2D k • p model, 
which not only captures correctly the evolution in band 
dispersion under a Lifshitz transition, but also describes 
the spin texture reasonabl}'^^^^. In this connection, it 
is natural to consider a minimal four-band model with 
Hamiltonian H{kx^ky) around X point on the [001] sur- 
face, which obeys the following three symmetries: mir- 
ror symmetry about the xz plane (M^z)^ mirror symme- 
try about the yz plane (Myz)^ and time-reversal sym- 
metry (© = TK, where K for denotes complex conju- 
gate). We then have under these symmetry operations, 

M^,H{k^,ky)M-^ = H{k^,-ky), My,H{k^,ky)M-^ = 

H{-k^, ky) and TH{k:,, ky)T-^ = i/*(-/c^, -ky). 

As shown in Fig. 2a, two distinct Dirac points occur 
at X, associated with energies Ej^ > Eq and E- < Eq. 
[Eq denotes the energy level at which the two cones in- 
tersect.] To account for the doublet state at each Dirac 
point, we choose the basis set to be the eigenvectors of 
Myz with eigenvalues ±i: {|i;Sn), | — z;Sn), | — i;Te), 
|i;Te)}, where the main atomic portion for each cone 
is indicated. In particular, when combined with the 
dominant orbitals for Sn (p^) and Te {px) atoms men- 
tioned earlier, this basis set also captures spin informa- 
tion, resulting in {|p;,,^;Sn), |p^,^;Sn), |p^,^;Te), 
|pa^,^;Te)}. Note that the quantization axis for spin is 



now along x, with Myz\pz^ -^ (^)) = ^^Pz^ -^ (^)) ^^^ 
Myz\px, -^ (^)) = Ti\Px^ -^ (^)). Defining 4x4 matri- 
ces, Sq,^ = 5q,(8)<J/3, with Pauli matrices s and a acting on 
spin and orbital spaces, respectively, the symmetry op- 
erator Myz then takes the form, iT^ss and the other two 
symmetry operators can be written as Mxz = —^^20 and 
T = —iTj20' After examining the 16 H matrices under all 
three symmetry operations, up to linear coupling in kx^ 
ky^ we obtain the following symmetry- allowed Hamilto- 
nian, 



H{k) = mHos + rn^^22 

+ kx{vix^20 + ^2x^02 + ^3x^23) 

+ ky{viyT^so^V2yT^ii^vsyT^ss). 



(1) 



In Fig. 2b, we plot the energy dispersions based on 
this effective Hamiltonian, which is clearly seen to cap- 
ture the key features of first-principles calculations (See 
Fig. 2g). The presence of two new Dirac points is due 
to the double degeneracy given by different Mxz eigen- 
values and is thus protected by the mirror symmetry of 
the system about the xz plane, in sharp contrast to the 
Dirac points at X, which are mainly protected by the 
time-reversal symmetry. As to the spin texture, note 
that (53) and (52) now represent the in-plane spin x and 
y components, respectively. Furthermore, one can prove 
that for any eigenstate with momentum k, (k|5i|k) = 
(out-of-the-plane component), as required by the com- 
bined TRS and the two mirror symmetries: (k|5i|k) = 
{k\MxzMyzesie-^M-,'M-^^\k) = -{k\si\k). It turns 
out that the resulting spin texture is qualitatively simi- 
lar to that shown in Fig. 3. Finally, we note that if one 



applies a unitary transformation, U ■ 
(1), the transformed H{k) becomes 
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where Vx± = vix ± V2,x and Vy± = Viy ± V2,y. This is in 
fact a more familiar Hamiltonian, representing two in- 
teracting Dirac cone s in t he Rashba form, k x s, with 
an identical chiralit}^^'^. In Table. I, we list the pa- 
rameter sets used in Eq. (2) by fitting experimentally 
obtained ARPES band dispersions in topologically non- 
trivial (Pb,Sn)Tei^ SnTe'*^^ and (Pb,Sn)S^^^^. 



IV. DISCUSSION AND CONCLUSION 

In general, topological crystalline insulators should 
harbor distinct classes with positive and negative mir- 
ror Chern numbers. The schematic diagrams of possible 
spin textures around the lower Dirac cones are given in 
Fig. 4. In Fig. 4a, a pair of Dirac cones with counter- 
clockwise spin texture on the horizontal mirror line gives 
the mirror Chern number n^ = — 2. The opposite case 
with a clockwise spin texture in Fig. 4b would give the 
mirror Chern number n^ = 2. By comparing the spin- 
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FIG. 4: (Color online) Schematic diagrams for spin tex- 
ture of the lower Dirac cones in distinct topological 
phases associated with various values of the mirror 
Chern number Um. a. ^m — —2; b. nm = 2; c. Hm — —1; 
d. nm = 1. Blue lines indicate the axis of mirror symmetry. 



texture obtained from our first-principles calculations in 
Fig. 3, we can conclude that the mirror Chern number 
of SnTe is —2. For Bi2Se3 system shown in Fig. 4c, on 
a mirror line, only a single Dirac cone appears at the 
center of the hexagonal Brillouin zone with the mirror 
Chern number n^ = — 1. In Fig. 4d, if the spin texture 
runs clockwisely, the system has the mirror Chern num- 
ber n^ = 1. Since there are two coaxial Dirac cones in 
SnTe, where the Sn-like and Te-like Dirac cones exhibit 
the same spin charality as Bi2Se3, it follows that the mir- 
ror Chern number of SnTe is twice that of Bi2Se3. 

In conclusion, we have delineated the charge density 
distributions and spin textures of the [001] surface states 
on the topological crystalline insulator SnTe via first- 
principles calculations. We show that the SnTe surface 
states can be pictured as two interacting coaxial Dirac 
cones which intersect each other to form the Dirac points 
along the mirror symmetry line f X. From an examina- 
tion of the charge distribution, we attribute the origin 



of the two coaxial Dirac cones to the two distinct atomic 



TABLE I: Parameter sets for the 2D k • p model Hamiltonian 
of Eq. (2) obtained by fitting the e xperi mental ARPES b and 
dispersions in (Pb,Sn)Tji^, SnT^^i"^ and (Pb,Sn)Se^324l 
{m, m , } are given in units of eV, and {vx±, V2x, Vy±, V2y} 
in units of eV-A. 
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species, Sn and Te, in the material. The crystal and time- 
reversal symmetries guarantee that the out-of-the-plane 
spin polarization is zero. The spin texture is dictated by 
the non-trivial Chern number n^ = — 2 in SnTe, which is 
different from the known Z2 topological insulators such 
as Bi2Se3. 
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